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Abstract 



We present a simple method to deal with caustics in the semiclassical 
approximation to the thermal density matrix of a particle moving on the line. 
For simplicity, only its diagonal elements are considered. The only ingredient 
we require is the knowledge of the extrema of the Euclidean action. The 
procedure makes use of complex trajectories, and is applied to the quartic 
double-well potential. 
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I. INTRODUCTION 



In the path integral formulation of quantum statistical mechanics, the thermal density 
matrix p(x,x') = (x\exp(—(3H)\x f ) of a system with Hamiltonian 



is given by 



p(x,x') 



*(W=* / s[z] 

[Vz{r)\ exp U 

z(o)=x' V n 



with 



S[z] 



dr 



- mz 2 + V(z) 
2 



(i) 



(2) 



(3) 



A semiclassical series for p(x, x') may be obtained from Eq. (|2]) through the method of 
steepest descent. The derivation depends solely on the knowledge of the paths that are 
minima of the Euclidean action S (the Euclidean nature of the path integral allows us to 
discard saddle-points). They act as backgrounds upon which a semiclassical propagator can 
be obtained exactly and then used to construct the series perturbatively. Its first term is 
given by 




(4) 



The sum runs over all minima x 3 c (t) of the action S[z] satisfying the boundary conditions 
z(0) = x' and z(/3h) = x, 



and Aj denotes the determinant of 



F[x; 



-m 



d 2 

— + V" 

dr 2 



(5) 



the operator of quadratic fluctuations around x{{j). (A derivation of this result will be 
sketched in Sec. |I[) 

In previous works @,|], we presented the explicit construction of the series for the di- 
agonal elements of the density matrix, p(x, x). For the sake of simplicity, we restricted our 
discussion to potentials of the single-well type. The more intricate case of multiple-wells - 
of which the quartic double well, with its many applications of practical importance 0, is a 
paradigm — was left aside, as it requires special treatment. Differently from single- wells, for 
multiple- wells the number iV of minima of S depends on x and (3 0. On the frontier sepa- 
rating regions in the (x, /3)-plane with different values of iV — a caustic — p sc (x, x) diverges 
due to the vanishing of the fluctuation determinant around the minimum that appears or 
disappears there. In the present context, this divergence is an artifact of the semiclassical 
approximation. Thus, a simple manner of eliminating it is certainly called for; this is the 
purpose of this paper. 
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As the caustic problem appears in other contexts in physics, it is instructive to briefly 
review how it comes about, and how it has been dealt with, for the sake of comparison. In 
optics, caustics occur whenever light rays coalesce. Thus, they separate regions of different 
number of extrema (the light rays) of the optical distance (the analogue of the action). 
In order to go beyond geometrical optics, one has to take into account fluctuations around 
these light rays. Just as in the present case, singularities emerge when we compute quadratic 
fluctuations on caustics. Ways to avoid this have been known for some time P-[T^] . Indeed, 
due to the traditional analogy between wave optics and quantum mechanics, the techniques 
involved are similar to the ones used in deriving connection formulae for WKB approxima- 
tions [EJ, and consist essentially in replacing one or more of the Fresnel integrals that arise 



in the stationary phase approximation with a so-called diffraction integral, whose form is 
specified by the classification of the caustic according to Catastrophe Theory [0,0 ; in the 
simplest case it is an Airy-type integral. A general procedure has also been developed to 
deal with caustics in the path integral formulation of quantum mechanics [|HJ. Although 
this general procedure could, in principle, be adapted to the case at hand, the nature of our 
problem allows for simplifications which warrant special treatment. 

In nonequilibrium quantum statistical mechanics, caustics are also known to occur in 
semiclassical descriptions of the decay of metastable states. The problem here is that of a 
particle in a potential which has a local minimum separated by a barrier from a region where 
it is unbounded below, in contact with a thermal reservoir. The phenomenon of caustics has 
been associated with a transition from the classical to the quantum regime of the decay rate 



T7H19| • General prescriptions for dealing with this phenomenon near the top of the barrier 



have been given in great detail [[I9|-|2T 



The case we shall analyze in this article differs from the one in the previous paragraph in 
the following aspects: (i) we discuss a problem in equilibrium quantum statistical mechan- 
ics; (ii) our analysis is global, in the sense that we compute the density matrix diagonal for 
every point on the real axis. It also differs from analyses carried out in optics and quantum 
mechanics because of the Euclidean nature of the path integral: only minima are to be con- 
sidered; saddle-points are discarded. (This is strictly true only in the "usual" semiclassical 
approximation. Our "improved" approximation makes use of some of the saddle-points.) In 
fact, in the specific example we analyze (the quartic double- well potential) only one new 
(local) minimum is introduced in function space after the various catastrophes that occur 
as we change the temperature; as it appears after the first catastrophe, this is the only one 
we have to consider. (Again, this is valid only for the "usual" semiclassical approximation; 
the improved one also requires an analysis of the second catastrophe.) This results in a 
prescription for dealing with caustics that is simple and direct. 

Previous works have studied the quartic double- well potential at finite temperature using 
semiclassical methods ||2^]23|| . Our work complements and extends those studies, by giving 



an explicit recipe for dealing with caustics. Variational methods have also proven extremely 
useful in this problem, and were quite successful in addressing applications in condensed 
matter physics [23]. Combinations of perturbation theory with variational techniques have 



also been recently used |25]. Our contribution to the semiclassical treatment opens the 
way for practical calculations, to be compared with perturbative, variational and numerical 
results. 

This article is organized as follows: before introducing the improved semiclassical ap- 
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proximation, which will remedy the problem of spurious divergences, Sec. [II] briefly reviews 
the derivation of the usual semiclassical approximation to the density matrix. The new 
method is, then, presented in two alternative ways: one has a better physical motivation, 
and clearly illustrates the essential ideas, but lacks effectiveness as a calculational tool; the 
other provides a general recipe to perform the calculations in a systematic way, resorting 
to the use of complex trajectories. Sec. |T| describes in detail the results obtained for the 
case of the quartic double- well potential by using the improved semiclassical approximation, 



which are compared to the usual approach. Sec. [TV] presents our conclusions and points out 
directions for future work. 



II. IMPROVING THE SEMICLASSICAL APPROXIMATION 

A. The "usual" semiclassical approximation 

In order to show how one can improve the semiclassical approximation so as to eliminate 
the unphysical divergences at the caustics, it is convenient to remember how the usual 
semiclassical approximation to a path integral like the one in Eq. (Q) is derived. Briefly, one 
has to: 

(A) Solve the Euler-Lagrange equation, m'z = V'(z), subject to the boundary conditions 
z(/3h) = z(0) = x, and determine, among the solutions, those which minimize (globally or 
locally) the action. For simplicity, we shall assume for the moment that there is only one 
such solution, which we denote by x c ; 

(B) Expand the action around x c : S[x c + rf\ = S[x c ] + S 2 + 5S, where 

1 rP h 

S 2 = - dr^F^^r), (6) 



2 Jo 



00 1 rPK 

5S = En drV^[x c (r)} V k (r); (7) 

k=3 Kl J0 



F is the operator defined in Eq. (H), and we are assuming that Viz) is an analytic function 
of z, so that all derivatives V^(z) exist; 

(C) Expand the fluctuations 77 (r) in terms of the orthonormal modes of the fluctuation 
operator F[x c ]: 

00 

= (8) 
3=0 

where Fipj{r) = Xj(fj{r), with (pj(0) = (pj(/3h) = 0; then 

1 00 

S 2 = ;EV]. (9) 
z i=o 



5S 



00 00 

EE- 

n=3 ii=0 



00 



}_ C {n) 

n=0 



(10) 
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where 

^.in = rdTV^[x c (T)} Vll (T)... Vln (r). (11) 
J 

The "usual" semiclassical approximation is obtained by neglecting 8S in the path integral 
on the r.h.s. of Eq. (^), which, upon the change of variables ry(r) — >■ {%}, becomes a product 
of Gaussian integrals: 

p(x,x)ne W (-^^f[l j , (12) 



Hence 

p(x,x)^exp^-^l^ A" 1 / 2 , (14) 

where A = E[^=o <\j = det F. (Explicit expressions for A [see Eq. ([|1]) below] were derived 
in Ref. j|, where it was also discussed how to systematically include corrections due to SS.) 
If there are N minima, one has to add together their contributions, thus obtaining Eq. (H). 



B. Taming the caustics 

When we cross a caustic, a classical trajectory x c (t) is created or annihilated. Precisely 
at this point, the lowest eigenvalue of F[x c ] vanishes, thus making the integral X blow up. 
This problem can be remedied by retaining fluctuations beyond quadratic in the subspace 
spanned by <po (the eigenmode of F associated with Ao), i.e., we replace X with 

~ [°° da ( V(oo)\ i /2 ( . 



where 

1 . , ™ 1 



V(a ) = -\ al + Y.-C£ ) ) ..o<- (16) 

A n=3 n - 

We take for M the smallest even integer such that C^} G is positive for all values of xq and 



(3; this suffices to make the integral in (ITS) finite even when Aq vanishes 



As a result, we obtain an improved approximation to the density matrix element (fj): 

p(x,x) » exp (-^A) A~ 1/2 T. (17) 



Here, x gm is the global minimum of S[x] 
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It is important to note that there is a one-to-one correspondence between the minima of 
S[x] and the minima of V(a ). Therefore, it is not necessary to explicitly add their contri- 
butions as in Eq. (f|), for they are already included in ([T7|) . Indeed, let a i(= 0), a 02 , ■ ■ ■ , a jv 
be the minima of V. If they are sufficiently far apart, one may compute T using the steepest 
descent method, obtaining 



N 




V(°oj) \ 



Substituting this into Eq. QT7| ) then yields 



(18) 



N 



S; 



p(x,x)^Eexp^-^)A~ i/i , (19) 



where Sj = S'fxgm] + V(a 0j ) and Aj = (V"(a j)/\o) A. Although expression ( |T9D is not 
identical to the "usual" semiclassical approximation [Eq. (f|)], the dominant term in both 
sums is the same, namely exp(— S[x gm ]/h) A -1 / 2 (recall that a 01 = 0, hence V(aoi) = and 
V"(aoi) = Ao); the other terms are exponentially supressed in the classical limit % — > 0. 

Another point that is important to mention is that one can, in principle, systemati- 
cally improve the "improved" semiclassical approximation, Eq. flTTD . To do this, one first 
decomposes the action into three pieces: S[x gm + rj] — S[x gm ] + Si + Su, where 

1 oo 

Si = V(oo) + ■= £ A i0 J (20) 



and 

5*11 = SS 



V(ao) - - A ao 



(21) 



with 5S defined by Eq. ( ^0|) . Applying this decomposition to Eq. (|) (with x' = x) then 
yields 

/ \ ( S[x gm ] \ [" d aj ( Si\~l( S u \ n 
p(x,x) = exp ^- / ]J -7= exp -- £ - -— . (22) 



This defines an "improved semiclassical series," the first term of which corresponds to Eq. 
(PI). Higher order terms can be readily computed, as they can be recast as sums of products 
of simple integrals. Compared with the "usual" semiclassical series 0, the series (^) has 
the disadvantage that integrals involving powers of ao must be computed numerically. On 
the other hand, those integrals are finite even at the caustics, so that the coefficients of the 
series (p2|) are well defined for any f3 and x . 

Although the procedure outlined in this subsection teaches us how to cross the caustics, 
it is not very convenient: in order to obtain the coefficients of V(ao) one has to find Ao 
and (fo{r). This, in general, is not an easy task, and makes the whole procedure very 
cumbersome. Instead, we shall present an alternative way of obtaining those coefficients, 
which is based on the one-to-one correspondence between the minima of S and V. 
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C. An alternative procedure 



Let us assume that M = 4 in Eq. (|i~6"D; this is the case for the quartic double- well 
potential, to be discussed in the next section. Then the "effective action" A(ao) = S[x gm ] + 
V(a ) for the "critical" mode <p is a fourth degree polynomial in a . 

Let us also assume for the moment that A(a ) has three extrema: a global minimum at 
ao = 0, a local maximum at u > 0, and a local minimum at v > u. This allows us to write 
A(a ) as 



A(a ) = S[x gm ] + a 



1 



UV <2q 



[U + V ) Oq + 



(23) 



(one can easily check that A'(0) = A'{u) = A'(v) = 0). 

We now have to relate a, u, and v to calculable quantities. We do this by imposing that 
A(v) = S[xi m ] and A(u) = S[x sp \, where xi m (r) and x sp (r) are the local minimum and the 
lowest saddle-point of S[x], respectively. This yields 



gmj 



(24) 



A(v)-A(0) e 3 (2-Q 
S[x sp ] - S[x gm ] A(u)-A(0) 2C-1 ' 

where ^ = v/u. It follows from the definition of x gm , x\ m and a; sp that the l.h.s. of Eq. (| 
is in the range [0, 1]. A plot of its r.h.s. shows that Eq. (|24l) possesses a unique real solution, 
lying in the interval 1 < £ < 2. 

Having determined £, we can now fix another combination of parameters, namely fi = 



au 



4. 



S[x sp ) - S[x gm ] = A{u) - A{<S) 



12 



We can then rewrite (f23|) as A(a ) = S[x gm ] + V3(a /u), where 



V 3 (z) 



2^ 



o^ 3 + ^ 4 



(25) 



(26) 



There still remains one parameter to be determined, namely u. Fortunately, we do not need 
it in order to compute T . Indeed, identifying Vz(ao/u) with V(ao) yields Ao = ^/u 2 , so 
that 



2nhu 2 



ddQ exp 



h 



(27) 



Changing the variable of integration to z = a^/u eliminates the unknown parameter u from 
the problem, leaving us with an expression for T which depends only on the calculable 
parameters £ and [i: 




(28) 



The case in which 5* has only one extremum can be dealt with similarly. Now ^4'( Q o) has 
one real root (a = 0), corresponding to the minimum x gm (r) of S[x], and a pair of complex 
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conjugate roots, w and w*, corresponding to the pair of complex conjugate trajectories x ct (r) 
and x* t (r). Correspondingly, we have A(a ) = S[x gm ] + Vi(a /\w\), where 



12 1 

- z 2 (cos 6) z 3 + - z 4 

2 3 V y > 4 



Vi(z) = X 

with x = ct\w\ 4 and = arg(u>). Identifying A(w) with S[x ct ] yields 



S[x ct ] - S[x 



X_ 
12 



2e 2i4> _ e a<j> 



(29) 



(30) 



from which we can obtain x an d </>• Finally, identifying Vi(ao/|w|) with V(ao) yields Ao 
x/|w| 2 , which leads to 



X 



2iTh\w\ 



dao exp 



Vi(a Q /\w\Y 



X 
2nh 



/ dz exp - 

J — oo \ 



Vi(^) 



(31) 



In the jargon of Ref. 0, the solution we have called x\ m {r) is a one-saddle: the operator 
of quadratic fluctuations around it, F[xi m ], has only one negative eigenvalue. At the second 
catastrophe, the one-saddle becomes a two-saddle (i.e., it becomes unstable in another di- 
rection in the functional space), and two one-saddles appear. Having lower action than the 
two-saddle, the one-saddles have a larger weight in the partition function. Hence, when the 
second caustic is crossed the role of "lowest saddle-point" in Eqs. (|24|) and (|25|) is transferred 
to one of the newly born one-saddles, namely the one with the lowest action. The transi- 
tion is smooth as all three trajectories coalesce and thus are identical to each other at the 
caustic. We also notice that, in spite of the infinite number of catastrophes, such a change 
of roles occurs only once, namely at the second catastrophe, since the minima of S and the 
one-saddles do not take part in the subsequent catastrophes. Indeed, as shown in Ref. 
(see also Sec. [Ill B 5| ), only (n — l)-saddles and n-saddles take part in the n-th catastrophe. 



III. APPLICATION: THE QUARTIC DOUBLE- WELL POTENTIAL 

A. Preliminaries 

Let us consider the quartic double-well potential, given by 

V(x) = ^ (x 2 - a 2 ) 2 (A>0). (32) 

In order to simplify notation, it is convenient to replace x and r by q = x/a and 9 = ujt, 
respectively, with uj = (Aa 2 /m) 1 / 2 . In the new variables, the equation of motion reads 
q = U'(q), where U(q) = \ (q 2 — l) 2 . Its first integral is 

1 -q 2 = U(q)-U(q t ), (33) 

where q t denotes the turning point (i.e., the point where q = 0). This can be further inte- 
grated to give us the relation between q t and the initial position q Q for a given "time of flight" 
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G = /3Hu. (As shown in Ref. |4[], that relation is all we need in order to compute the "usual" 
semiclassical approximation to p(x,x), Eq. (£|). The same is true for the computation of the 
improved approximation, Eq. (|TT|), using the procedure outlined in Sec. [11 C .) Assuming for 
the moment that < q$ < q t < 1, we can write 



6 
~2 



dq 



w J2[U(q)-U(q t )\ 



(34) 



Inserting the explicit form of U(q) and changing the integration variable to z = q/qt, Eq. 
(1341) becomes 



u 



dz 



qo / qt ^{\- z 2 )(l-k 2 z 2 ) 



(35) 



where 







u = 



k 2 = 



Qt 



2~q( 



2 ' 



(36) 



Performing the integration (formula 130.13 of Ref. pTfl ) and solving for q finally yields 

q = q t cd(u,k), (37) 

where cd is one of the Jacobian elliptic functions. 

The action can be written as S[x] = (h/g) I[q], where g = h\/m 2 uj 3 and 



dB 



1 -2 

- Q 

2 H 



U(q) 



Using Eq. (p3l), we may rewrite I[q c ] as 



I[q D ] = 6 U(q t ) + 2 r dq p [U(q) - U(q t )}. 



(38) 



(39) 



The integration can be done with the help of formula 219.11 of Ref. [27|. After a few 
algebraic manipulations one arrives at 



I[Qc 



QU(q t 
2 
3 



y2q 2 (q 2 -q 2 )(2-q 2 -q 2 ) 
-\ p(2-q 2 ){(l - q 2 t ) [K(k) - F( V , k)] - E(k) + E(^, k)} 



(40) 



where K, F and E are elliptic integrals 0,^1 anc ^ V 9 = arcsin(go/?t)- 

Equations (^7| ) and fl40|) have been derived under the assumption that go and q t are 
real and satisfy < qo < q t < 1. However, since the elliptic functions and integrals are 
meromorphic functions of their arguments, we can now abandon that assumption and treat 
go and q t as complex variables. (Note, however, that I[q c ] is a multivalued function of q t and 
so one must be a bit careful when computing it. For instance, the first square root in Eq. 
(f!0|) acquires a minus sign if — 1 < q t < —q < 0.) 
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Finally, the determinant of the fluctuation operator is given by [[| 



A = 4„sgn( 90 - 9t ) ^ ^—j^. (41) 

We now have all the ingredients to compute the semiclassical approximation to p(x, x) 
- both the "usual" and the "improved" one. Indeed, both the action and the determinant 
of fluctuations can be expressed solely in terms of q t . Therefore, as anticipated, this is the 
only information we need from the classical trajectories. 



B. Singularities and their removal 

As we have already said, the "usual" semiclassical approximation to p(qo, go) diverges at 
a caustic because of the vanishing of the determinant of fluctuations A around the minimum 
of S which appears or disappears there. According to Eq. (f|l]), there are two ways A may 
vanish: (i) when (dq /dq t )s = 0; (ii) when U(q ) = U(q t ). A qualitative analysis of the 
equation of motion shows that, at the boundary between the N = 1 and the N = 2 regions 
in the (g , 0)-plane, A vanishes according to the first alternative JFj. Solving the equation 
(dqo/dq t )e = for q t and inserting the result g~t(0) into Eq. (|37|) , one obtains the lower 
curve depicted in Fig. ^ — the caustic. 

In what follows we shall examine the behavior of the "usual" semiclassical approxima- 
tion across the caustic, and compare it with the improved approximation. (All numerical 
calculations were performed using maple.) 



1. q = 0, 6 « tt 

When go — and < tt, the only real solution to Eq. (^) is q t = 0. It then follows from 
Eq. ( |39|) that I[q c ] = 0/4. In order to compute A we need qo(q t , O) for small q t . Using Eqs. 
(§§) and Q37p we find q w g t cd(6/2,0) = q t cos(6/2); Eq. (^j) then yields A = 2ng sin 6 



in the limit q t — > 0. Therefore, the "usual" semiclassical approximation to p(0, 0) gives 

p sc (0,0) « (2ng sin6)- 1/2 exp ( J (6 < vr). (42) 



It diverges like (tt — 0) _1//2 as ^ tt~ . 

While for 9 < tt there is only one real solution to the equation go(<Zt> ©) = 0, for > tt 
there are three: q t = 0, corresponding to the trajectory q c (0) = (which is now a 1-saddle), 
plus a pair of solutions located symmetrically with respect to the origin, corresponding to 
a pair of degenerate minima of the action (see Fig. |^). The latter can be traced back to a 
pair of purely imaginary trajectories for < tt. Indeed, making q t = i£ in Eqs. (|3"6f) - 



and using the identity cd(u, ik) = cn(w a/1 + k 2 , k/ y/1 + k 2 ) p7| , we obtain 

qo = itcn (| JlTl 2 , . } ) . (43) 
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The r.h.s. of the above equation has an infinite number of zeros besides the one at £ = 
(see Fig. |3[). As approaches ir from below, the zeros approach the origin and two of them 
eventually coalesce there when = 7r, reappearing as a pair of real solutions to the equation 
9o(9t,Q) = for > 7T. 

In Fig. | we show both the usual [Eq. (|])] and the improved [Eq. (|T7|)] semiclassical 
approximation to p(0, 0) for « 7T. 

jg. G = vr, g ~ 

When = 7r, the approximation g ~ <7t cos(0/2) is not enough for our purposes. Going 
to the next nontrivial order in the Taylor expansion of qo(qt, vr) one obtains go ~ Qt as It ~* ®- 
It then follows from Eqs. (f|), ([4(]) and fl4T| ) that the usual semiclassical approximation to 
p(qo, qo) behaves, for = tt, as 

p sc (%, qo) 9- 1/2 \qo\- 1/3 exp (-^) . (44) 

Two aspects of this result are worth of mention: (i) the singularity at q Q = is integrable, 
hence the semiclassical partition function is well defined; (ii) because of the exponential 
factor, if g -C 1 one has to be very close to the origin to "notice" the singularity: for 
p sc (qo, qo) to be of order unity or greater, q must satisfy |g | ~ <? 2 ex P(~ 37r/4^) . 

Fig. ^ shows both the usual and the improved semiclassical approximation to p(qo, qo) for 
= 7T. In order to make visible the singular behavior of the former, we have taken g = 0.3. 
One can notice that far from the caustic (i.e., for q > 0.2) the two curves are similar, but 
differ by approximately 10%. This is due to the relatively large value of g. The difference 
becomes smaller as g is made smaller, and vanishes in the classical limit g — > 0. 

s. e > vr 

Expanding q (q t , 0) around q t (Q), we obtain q -q~o ~ (<?t-<?t) 2 , so that dq /dq t ~ (g t -?t) 
near the caustic. The other terms in Eq. (|41| ) remain finite on it, so that we finally obtain 

A- 1/2 ~|g t -g t |- 1/2 ~|g -g r 1/4 . (45) 

Fig. |] depicts both the usual and the improved semiclassical approximation to p(qo, qo) 
for = 5.0. Again we had to use a relatively large value of g in order to magnify the 
"critical" region where the usual semiclassical result diverges. Note that the divergence 
occurs only at the two minima side of the caustic (see Fig. 0), as it is associated with the 
coalescence of the local minimum with a saddle-point of the action; the contribution of the 
global minimum remains finite at the caustic. 

4. G > 2vr 

As discussed in Ref. [0], another catastrophe is present if > 2tt. This time, A vanishes 
when the classical trajectory is such that U(q ) = U(q t ), or, since the potential is symmetric 
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about the origin, q t = — g . This catastrophe is associated with the appearance of periodic 
classical trajectories, and the condition q t = — g determines the amplitude A(Q) of these 
trajectories [^4(0) is the positive solution to equation qo(qt, ©) = —qt, where go(gt,0) is 
defined by Eqs. (|36| ) and (|37|)1. It is not difficult to see why a catastrophe occurs when that 
condition is satisfied: if |go| < ^4(0), there are two periodic trajectories satisfying q c (0) = qo, 
related by time reversal, i.e., qj?'(0) = q^i® — 9). If, on the other hand, |go| > A(Q), no 
such trajectories exist. |go| = A(Q) thus marks the boundary between regions with zero and 
regions with two periodic trajectories. This boundary is depicted in Fig. [I] (upper curve). 

As discussed at the end of Sec. [11 Q the procedure for dealing with the caustic outlined 
in that section is not affected by the appearence of a new catastrophe. What changes as the 
second catastrophe is crossed is the identity of the "lower saddle-point": for \q \ > A(Q), it 
is to be found among the solutions of Eq. (|3~7D; for \q \ < A(Q), it is given by any one of 
the two periodic trajectories satisfying g c (0) = g (since they have the same action). As a 
matter of fact, since all periodic trajectories with the same amplitude (and the same period) 
have the same action, we may pick the one that satisfies the condition q t = —qo, for then 
we can use Eq. (pEO|) to compute its action. 



5. 6 -» oo 

In the high-temperature limit, G — > 0, thermal wavelengths are very small, and the 
classical limit sets in, as quantum fluctuations are suppressed. That is the regime where our 
improved semiclassical approximation should work best, as illustrated in Ref. since it 
incorporates quantum fluctuations in a controlled manner as we lower the temperature, and 
profits from the simplification of having to deal with only one or two minima, as already 
emphasized. Nevertheless, even in the opposite limit, — > oo, our improved semiclassical 
method can be used to reproduce zero temperature results. The secret is to recognize that 
the various saddle-points which were discarded for finite do play a role in such a limit; in 
fact, taking them into account is equivalent to using a dilute-gas approximation, as we will 
qualitatively argue. First, however, let us review how the saddle-points emerge. 

For a fixed g in the interval (—1, +1), new saddle-points will appear as we increase 9, 
following a pattern outlined in Ref. @. As a result, the strip of the (go, ©)-plane defined 
by — 1 < g < +1 and > may be divided into regions wherein each (g , 0) point gives 
rise to 2n + 1 solutions, n = 0, 1, 2, . . ., as shown in Fig. [I]. The regions are separated by 
caustics where instabilities develop: for n even, as we cross the caustic between the regions 
with 2n + 1 and 2n + 3 solutions, a n-saddle, i.e., a solution with n negative eigenvalues, 
and a (n + l)-saddle appear; for odd values of n, the n-saddle in the region with 2n + 1 
solutions becomes unstable — it is replaced, in the region with 2n + 3 solutions, by two 
n-saddles which are periodic, and by a (n + l)-saddle. Thus, for a given n, we end up with a 
(n + l)-saddle, (n — 1) pairs of n-saddles with n < n, and a pair of minima (two 0-saddles). 

As — > oo, the two minima will correspond to solutions q c (9) that spend most of their 
euclidean time 9 near either go = —1 or go = +1. They both have a single turning point, 
and only the local minimum will cross the origin q = 0: first, at a small value of 9] and 
upon returning, at a large value 9 ~ 0. The 1-saddles, which are periodic, have two turning 
points, will also cross the origin twice, but one of the crossings will occur at a value of 9 
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near 0/2. Generalizing this qualitative analysis, we may conclude that n-saddles will have 
(n+ 1) turning points and n inner (away from 8 = and 9 = 0) crossings of the origin. Each 
of those inner crossings is equivalent to a kink or an antikink, so that a generic n-saddle will 
not differ much from a solution built out of a superposition of kinks and antikinks in that 
limit. Typically, the euclidean time width of such kinks should be much smaller than their 
separations, as the solutions tend to spend most of their time near go = ±1- 

Varying the euclidean times where those inner crossings occur should not alter signifi- 
cantly the action of the n-saddle in the — > oo limit, an indication of the existence of a 
flat direction in functional space corresponding to that variation. As flat directions are as- 
sociated to near-zero eigenvalues, we claim that the negative eigenvalues which characterize 
the n-saddle will approach zero from below as — > oo, and that the euclidean times where 
the crossings occur should be treated as collective coordinates [p8|l , just as the positions of 
kinks and antikinks in the dilute-gas approximation. The contributions of the various kinks 
and antikinks can be dealt with in the usual manner — they add up to an exponential, and 
reproduce the standard result for the splitting between ground and first-excited state |29 



(See, however, Ref. [[HJ for a treatment of the low-temperature limit within the functional 
formalism that does not appeal to the dilute-gas approximation.) 



IV. CONCLUSION 

Semiclassical methods are a powerful nonperturbative tool, for both equilibrium and 
nonequilibrium systems. This article, together with Refs. [HJjl[7j, represents a further step 
towards a systematic semiclassical treatment of quantum statistical mechanics. 

In the present work, we developed a simple procedure to derive the lowest order semiclas- 
sical approximation for the case of multiple-well potentials in equilibrium quantum statistical 
mechanics. In order to adequately incorporate new extrema, we kept fluctuations beyond 
the quadratic level along the "unstable" direction in functional space, and relied on our 
knowledge of the type of catastrophe involved as we cross a caustic to eliminate spurious 
singularities in the semiclassical approximation, obtaining sensible results for the density 
matrix elements for any temperature. This was exemplified by the analysis of the quartic 
double-well potential. 

Our results can possibly be extended to nonequilibrium systems, such as those where 
a time-dependent potential is coupled to a heat bath, in order to better understand tran- 
sient regimes. Although the physics of nonequilibrium quantum statistical mechanics has 
been considered in detail in the context of semiclassical calculations of the decay rates of 
metastable systems JTT|-|2"T|] , a thorough analysis of the various transient regimes, and of 
the interplay of their corresponding time scales, is still needed. Here, however, we will no 
longer profit from the drastic reduction in the number of extrema that occurs in equilibrium 
situations, as time evolution forces us to deal with saddle points and maxima, as well. The 
simplified methods presented in this paper will still be useful to describe the asymptotic 
imaginary time evolution corresponding to equilibrium, but not the real time evolution, 
which requires the traditional quantum mechanical treament. 

As for possible extensions to field theories, the methods developed in ||,[|] should be 
applicable to the evaluation of the effective potential in the presence of non-trivial back- 
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grounds (defects), as long as they depend on only one coordinate. This can be of use in a 
wealth of possible applications, and should help in the study of phase transitions and critical 
phenomena where such defects play a role. Cases such as the ones explored here and in J?|, 
which involve several extrema, still lack a field theoretic treatment. 
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FIGURES 




FIG. 1. The lower curve in this figure depicts the caustic for the quartic double- well potential. 
Below it the action has only one minimum; above it, the action has two minima. The cusp is located 
at the point (qo, 0) = (0, ir). A second catastrophe occurs at the curve in the middle: below it (but 
above the caustic) the action has a one-saddle (in addition to the minima); as the curve is crossed, 
the one-saddle splits into a two-saddle and a pair of one-saddles, the latter corresponding to a pair 
of periodic trajectories. The minimum of the middle curve is located at (qo,&) = (0, 2ir). Upon 
crossing the upper curve, whose minimum is located at (qo,Q) = (0, 3ir), the number of classical 
trajectories increases by two: the two-saddle splits into a three-saddle and a pair of two-saddles. 
Numbered arrows correspond to the first three subsections of Sec. QIB[ (This figure corrects Fig. 
5 of Ref. 0.) 
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FIG. 2. q (q t ,Q) [Eq. ©] for 6 = 2.0 (solid line), G = n (short-dashed line), and Q = 4.5 
(long-dashed line). For G < tt this function is one-to-one. For Q > tt and \qo\ sufficiently small 
there are three (or more) real values of q t corresponding to a given qo. 




FIG. 3. -%(£,9) [Eq. ©] for G = 2.0 (solid line), G = 2.5 (short-dashed line), and Q = vr 
(long-dashed line). 
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FIG. 4. p(0, 0) vs. for g = 0.3. Usual (dashed line) and improved (solid line) semiclassical 
approximation. 
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FIG. 5. p(qo,qo) vs. qo for = n. Usual (dashed line) and improved (solid line) semiclassical 
approximation. 
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FIG. 6. p(qo,Qo) vs. qo for = 5.0 and g = 0.3. Usual (dashed line) and improved (solid line) 
semiclassical approximation. 
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FIG. 7. q t vs. qo for = 5.0. If |go| < q* = 0.3332..., there are three real solutions to 
Eq. ( |37| ) (solid lines); the upper curve corresponds to the global minimum, the lower one to the 
local minimum, and the one in the middle to a one-saddle. At the caustic the last two coalesce, 
and reappear at the other side of the caustic (i.e., |go| > <7*) as a pair of complex conjugate 
solutions. Their real and imaginary parts are represented by the long-dashed and short-dashed 
lines, respectively. 
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